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Abstract Bat coronaviruses (CoV) are putative precur- 
sors of the severe acute respiratory syndrome (SARS) CoV 
and other CoV that crossed the species barrier from zoo- 
notic reservoirs into the human population. To determine 
the presence and distribution of CoV in Iberian bats, 576 
individuals of 26 different bat species were captured in 13 
locations in Spain. We report for the first time the presence 
of 14 coronaviruses in 9 Iberian bat species. Phylogenetic 
analysis of a conserved CoV genome region (RdRp gene) 
shows a wide diversity and distribution of alpha and 
betacoronavirus in Spain. Interestingly, although some of 
these viruses are related to other European BatCoV, or to 
Asian CoV, some of the viruses found in Spain cluster in 
new groups of « and f CoV. 


GenBank accession numbers: Partial sequences of RdRp gene of bat 
coronaviruses obtained in this study were deposited in GenBank 
under the accession numbers HQ184049 to HQ184062. 
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The emergence of infectious diseases is a major threat to 
global public health in this century (WHO, World Health 
Report —_http://www.who.int/whr/previous/en/index.html) 
and many of these new infectious human diseases are 
caused by viruses emerging from wildlife. In the last 
50 years, more than 30 new infectious human diseases 
have been identified (WHO, World Health Report 
http://www.who.int/whr/previous/en/index.html) including 
the Severe Acute Respiratory Syndrome (SARS). The 
aetiological agent of this disease was identified as a 
previously unknown coronavirus (SARS-CoV) [1] and 
BatCoV are putative precursors of SARS-CoV [2]. The 
outbreak of SARS-CoV and subsequent identification of 
two additional human coronaviruses (HCoV-NL63[3] and 
HCoV-HKU1 [4]) has drawn human and animal health 
attention to Coronavirinae subfamily, that includes three 
genera, Alphacoronavirus (« CoV), Betacoronavirus 
(6 CoV) and Gammacoronavirus (y CoV), replacing the 
classical groups 1, 2 and 3 [5] (http://talk.ictvonline.org/ 
media/g/vertebrate-2008/default.aspx). 

The relevance and possible re-emergence of the pan- 
demic SARS-CoV and other emerging viruses of zoonotic 
origin have activated surveillance systems of hazard agents 
in wild animals, including bats. As a result of these studies, 
bats have been described as putative reservoirs for some 
emerging viruses affecting humans [6]. BatCoV are puta- 
tive precursors of CoV affecting humans and mammals [7, 
8], including SARS-CoV [2] and other CoV that crossed 
the species barrier from zoonotic reservoirs into the human 
population [9]. In fact, association of some of these CoV to 
certain bat species has been suggested [10, 11], reinforcing 
the notion that there may be a relationship between some 
BatCoV and their hosts. Nowadays the presence of CoV 
has been shown in bats in China [12, 13], North and South 
America [14-17], Africa [18] and a number of regions in 
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Europe [11, 19-21] but not the Iberian Peninsula (Spain 
and Portugal), which is a bridge for European and African 
bat populations [22]. Thus, surveillance of wildlife reser- 
voirs of putative zoonotic CoV is necessary, not only to 
unveil the ecology of these viruses, but also to permit early 
detection of viruses that might pose a threat to human 
health. 

To determine the presence and distribution of putative 
zoonotic CoV in Iberian bats, 576 individuals from 26 bat 
species were captured and sampled in 13 different locations 
throughout Spain during 2004-2007 (Fig. 1). These sam- 
ples were collected in the context of bat rhabdoviruses and 
lyssaviruses Surveillance Program in Spain. Most of the 
sampled bat species are also distributed across Europe, but 
Eptesicus isabellinus is a meridional serotine bat restricted 
to North Africa and the Iberian Peninsula [23], and Myotis 
escalerai is endemic in the Iberian Peninsula. Bats were 
caught with mist nets mainly as they left diurnal roost and 
by hand with polyethylene butterfly nets within roosts. 
Oro-pharyngeal swabs (n=390) between 2004 and 2007 as 
well as faecal samples from individual bats (n=216) in 
2007 were taken before bats were released (Table 1). 

Oro-pharyngeal swabs collected between 2004 and 2006 
were preserved in | ml of lysis buffer (4 M GuSCN (Sigma), 
OS% N-lauroyl Sarcosine (Sigma), 1 mM dithiothreitol 
(DTT, Sigma), 25 mM Sodium Citrate and 20 pg/tube 
Glycogen (Boehringer Mannheim). Oro-pharyngeal swabs 
and faeces collected in 2007 were preserved in both | ml of 
lysis buffer and 1 ml of viral transport medium (VTM) 
(Eagle’s minimal essential medium (EMEM) supplemented 
with 10 Ul/ml of penicillin, 10 g/ml of streptomycin, 
160 jig/ml of gentamicin, 50 UI/ml of mycostatin and 1% of 


Fig. 1 Geographical location of bat capture sites in Spain. 1: A 
Corufia, 2: Lugo, 3: La Rioja, 4: Gerona, 5: Castellon, 6: Valencia, 7: 
Alicante, 8: Malaga, 9: Cadiz, 10: Sevilla, 11: Huelva, 12: Caceres, 
13: Menorca. The stars show the locations where positive samples 
were found 
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bovine serum albumin). All samples were frozen at —80°C 
before sending them to the Rabies Reference Laboratory at 
the Centro Nacional de Microbiologia, ISCIII in Madrid. 
Faecal samples were clarified by centrifugation. Total 
nucleic acid was extracted from a 200 jl aliquot of each 
specimen for PCR assays and the rest were stored to —80°C 
in two different aliquots. Final pellets were always resus- 
pended to 55 ul of water. 

A pan-coronavirus nested PCR was designed in the 
RdRp gene. A total of 5 jl of extracted RNA was added to 
45 ul of reaction mixture of OneStep RT-PCR kit 
(QIAGEN, Valencia, CA, USA) containing 200 uM 
dNTPs and 60 pmol of generic CoV-specific degenerated 
primers (forward 5'-CARATGAATYTIAARTAYGC-3' 
and reverse 5'-TGYTGWGARCAAAAYTCRTG-3’) and 
following manufacturer indications. Amplifications were 
carried out into thin-walled reaction tubes (Sorenson, 
BioScience, UT) in a PTC-200 (Peltier Thermal Cycler, MJ 
Research, Watertown, MA). Nested PCR amplifications 
were performed using 2 wl of first amplification product 
and 48 ul of reaction mixture containing 60 mM Tris-HCl 
(pH 8.5), 15 mM (NH4)2SO,, 200 uM dNTPs (Amersham 
Pharmacia Biotech, Piscataway, NJ), 3mM MgCl, 
35 pmol of generic CoV-specific degenerated primers 
(forward 5'-ATGGGWTGGGAYTAYCCIAARTG-3’and 
reverse 5/-ACRTTRTTYTGRWARTA-3’) and 1.25 U 
AmpliTaq DNA Polymerase (Perkin-Elmer Cetus, Nor- 
walk). Amplification product size of 512 nt was visualized 
by agarose gel electrophoresis and sequenced directly in 
both directions using an automated ABI PRISM 377 model 
sequencer. For phylogeny reconstruction, 
sequences were aligned together with others obtained from 
public genomic databases using the program CLUSTAL X 
(version 1.83) (Table 2). A Bayesian phylogenetic infer- 
ence was obtained using Mr Bayes version 3.1 [24] with 
random starting trees without constraints. For the analyses 
GTR substitution model, gamma estimation and two 
simultaneous runs of 10” generations were done, each with 
four Markov chains, and the trees were sampled every 100 
generations. Amino acid identity was calculated with 
MEGA 4 using the pairwise deletion option. The alignment 
comprised the same 396 bp of the RdRp gene used for the 
phylogenetic reconstruction. 

A total of 26 out of the 30 known bat species known for 
Iberian Peninsula were screened for CoV and 14 samples 
taken from 9 bat species, all included in the family 
Vespertilionidae, were positive for CoV RNA (Tables | 
and 2). Twelve of them were found within faecal samples 
(5.5%) of 7 different bat species in 6 locations and two 
were obtained in oral samples (0.5%) of 2 other different 
species in the same location (Table 1). In view of these 
results, it is not surprising that the presence of CoV RNA is 
significantly more frequent in faeces than in oral cavity 
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Table 1 Results of detection of CoV RNA in faecal or oral samples of bats collected in Spain 

Bat Species Faecal samples Oral samples Location Genus 
positive/no. tested positive/no. tested 

Barbastella barbastellus 0/4 0/2 3,4 

Eptesicus isabellinus 1/8 NA 10° B 

Eptesicus serotinus 0/7 NA 1, 12 

Hypsugo savii 2/26 0/10 9,49, 107" ot, 

Miniopterus schreibersii 0/2 1/71 2, 5, 6°, 7, 8, 9, 12, 13 a] 

Myotis alcathoe 0/1 NA 3 

Myotis bechsteinii 0/2 0/3 4,8 

Myotis blythii NA 1/11 6°, 7, 12 a 

Myotis capaccinii NA 0/14 5, 6, 13 

Myotis daubentonii 1/39 0/52 2,.3;.87) 11, 12 oO 

Myotis emarginatus NA 0/2 8, 12 

Myotis escalerai NA 0/15 11,7 

Myotis myotis 1/1 0/17 6, 7, 8°, 12 a 

Myotis mystacinus 0/5 NA 2,3 

Myotis nattereri 0/3 0/3 4 

Nyctalus lasiopterus 5/37 0/137 3, 32> 9, 10, 11 at 

Nyctalus leisleri 0/23 O/11 2, 3, 4, 8 

Pipistrellus kuhlii 1/4 0/6 8, 12° a 

Pipistrellus pipistrellus 0/3 0/1 3, 12 

Pipistrellus pygmaeus NA 0/1 12 

Pipistrellus sp. 1/29 0/5 1, 2, 4, 8, 122° oO 

Plecotus auritus 0/7 NA 1, 3,4 

Plecotus austriacus 0/7 0/10 4,11, 12 

Rhinolophus euryale NA 0/13 6, 7, 12 

Rhinolophus ferrumequinum 0/3 0/5 4, 8, 12 

Rhinolophus hipposideros 0/4 NA 4,12 

Rhinolophus mehelyi NA 0/1 12 

Total 12/216 2/390 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 a,p 


NA no samples available 


* Locations where positives samples were found 


> These samples were collected in different localities than other positive samples with the same number location 


(p<0.0001, FISHER EXACT TEST). It is of interest that 
none of the viruses has been found in oropharingeal and 
faecal samples of the same individual, when both samples 
were available. This fact may indicate either that the 
infection was at different stage in the different individuals 
at the time of sampling or that replication of virus may take 
place independently in the intestinal and respiratory tracts 
[25]. Most of the CoV RNA sequences found in faecal 
samples (83%) correspond to « CoV, the remaining two 
belonging to f CoV. All CoV RNAs from oral samples 
(100%) were found to contain viral RNA sequences cor- 
responding to « CoV. In agreement with all previous 
studies [7], none of the coronavirus detected in Spanish 
bats belong to group y. 

The phylogenetic analysis of Spanish BatCoV was 
performed using 396nt out of the 512 nt RT-PCR amplified 


fragments. 116nt fragment information was lost to allowed 
us include more sequences from other European countries 
and other continents deposited in GenBanK to perform a 
meaningful analysis. This small part of the RdRp gene has 
been previously used, and sufficiently represents the full 
gene information, for phylogenetic analysis of BatCoV [11, 
14, 20] 

The phylogenetic reconstruction showed 6 different 
lineages of Spanish BatCoV (Fig. 2). BatCoV A and B 
were closely related to other « BatCoV found in China 
[12], although they appeared to display certain genetic 
differentiation (Fig. 2). Myotis daubentonii-associated 
CoV H, and Pipistrellus-associated CoV K, clustered 
respectively with lineages 4 and 3 of « CoV previously 
described in Germany and are hosted by the same bat 
species or genera [11] (Fig. 2). BatCoV L was closely 
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Table 2 Data of interest related to the 91 coronavirus sequences used for the generation of the phylogenetic tree 
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Access no 


DQ249221 
DQ249219 
DQ249218 
DQ648809 
DQ648807 
DQ249217 
DQ648819 
DQ249215 
DQ249214 
DQ074652 
DQ249216 
DQ648803 
HQ184059 
HQ184062 
GQ404795 
GQ404796 
GQ404797 
DQ022305 
NC_009696 
NC_004718 
NC_009021 
NC_006577 
NC_006852 
NC_007732 
NC_005147 
EF544563 
EF544565 
HQ184049 
HQ184050 
DQ648838 
DQ648855 
DQ648854 
NC_003436 
EU375862 
EU375859 
EU375858 
EU375855 
EU375863 
EU375861 
EU375856 
EU375854 
EU375857 
EU375865 
EU375853 
EU375860 
EU375869 
EU375864 
EU375870 


Host species 


Bat 

Bat 

Pipistrellus sp. 

Bat 

Bat 

Pipistrellus sp. 

Bat 

Tylonycteris 

Tylonycteris 

Tylonycteris 

Tylonycteris 

Bat 

Hypsugo savii 

Eptesicus isabellinus 
Rhinolophus hipposideros 
Rhinolophus hipposideros 
Rhinolophus hipposideros 
Rhinolophus sinicus 
Rhinolophus macrotis 
Human 

Rousettus leschenaulti 
Human 

Mouse 

Pig 

Human 

Myotis occultus 

Myotis occultus 
Miniopterus schreibersii 
Myotis blythii 

Bat 


Rhinolophus ferrumequinum 


Rhinolophus ferrumequinum 


Pig 

Myotis dasycneme 
Myotis dasycneme 
Myotis dasycneme 
Myotis dasycneme 
Myotis dasycneme 
Myotis dasycneme 
Myotis dasycneme 
Myotis dasycneme 
Myotis dasycneme 
Myotis bechsteinii 
Myotis bechsteinii 
Myotis bechsteinii 
Pipistrellus nathusti 
Pipistrellus nathusti 


Pipistrellus pygmaeus 


Country 


China 
China 
China 
China 
China 
China 
China 
China 
China 
China 
China 
China 
Spain 
Spain 
Slovenia 
Slovenia 
Slovenia 
China 


China 


USA 

USA 

Spain 
Spain 
China 
China 
China 


Germany 
Germany 
Germany 
Germany 
Germany 
Germany 
Germany 
Germany 
Germany 
Germany 
Germany 
Germany 
Germany 
Germany 


Germany 


Genus 


RF RFR RBRRRKCRKCRTRTBDBDpWp pvp vv vo vvvvwvvwvovwvvwpvvwvovoovovwlvvvov0o»l0d 


Rg R8 RR QB RRR 


Rg RR Rk QB RRR 


Cluster 


Nn 


HKU 
HKU 
HKUS 


Nn 


HKUS 


HKU 
HKU 
HKU 
HKU 


i a 


HKU9 
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Table 2 continued 


Access no Host species Country Genus Cluster 
EU375868 Pipistrellus pygmaeus Germany o 

EU375867 Pipistrellus pygmaeus Germany a 

HQ184060 Pipistrellus sp. Spain w% 

DQ648822 Bat China a 

DQ648821 Bat China a 

DQ648824 Bat China a 

DQ648823 Bat China a 

EU375875 Myotis daubentonii Germany of 

EU375873 Myotis daubentonii Germany a 

EU375874 Myotis daubentonii Germany a 

EU375872 Myotis daubentonii Germany a 

EU375866 Myotis daubentonii Germany a 

HQ184056 Myotis daubentonii Spain oA 

EU375871 Myotis daubentonii Germany a 

DQ648833 Myotis ricketti China oA 

DQ249224 Myotis ricketti China a HKU6 
DQ648837 Myotis ricketti China oA 

DQ249235 Rhinolophus sinicus China a HKU2 
DQ249213 Rhinolophus sinicus China oA HKU2 
DQ648840 Bat China a 

DQ249228 Miniopterus sp. China a HKU8 
EU834954 Miniopterus australis Australia o 

EU834952 Miniopterus australis Australia a 

EU834955 Miniopterus schreibersii Australia ox 

EU834953 Rhinolophus megaphyllus Australia a 

DQ648835 Miniopterus schreibersii China oA 

DQ648796 Bat China 4 

DQ648797 Bat China a 

DQ249226 Miniopterus magnater China oA HKU7 
HQ184061 Hypsugo savii Spain a 

HQ184051 Nyctalus lasiopterus Spain oA 

HQ184054 Nyctalus lasiopterus Spain a 

HQ184053 Nyctalus lasiopterus Spain oA 

HQ184052 Nyctalus lasiopterus Spain a 

HQ184055 Nyctalus lasiopterus Spain oA 

HQ184057 Myotis myotis Spain oA 

HQ184058 Pipistrellus kuhlii Spain a 

EU834951 Myotis macropus Australia a 

NC_002645 Human a 

NC_005831 Human of 

AY994055 Cat OL 

NC_002306 Pig o 

NC_001451 Chicken y 
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0,64 


BICoV-HKUS (DQ249221) 
|BtCoV-HKUS (DQ249219) 
| 'BtCoV-HKUS (DQ249218) 


0,94, _ ‘Brcov (00648807) 
~~ | “BtCoV-HKUS (DQ249217) 
BtCoV (DQ648819) 
= BtCoV-HKU4 (DQ249215) 
Jexcov-HKU4 (00249214) CMO58 
|'BtCoV-HKU4 (DQ249216) 
“BtCoV (DQ648803) 
0,52 BtCoVd (HQ184959) 
‘BICoV M (HQ184962) + 83,1 
0.91{BtCov (Ga404795) 
=" BICoV (GQ404796) 
4|  ‘BtCov(Ga404797) 
—— 0,99 —SARS-like BtCoV (DQ022305) 
0,87 L SARS-ike BtCoV (NC_009696) 
= SARS-CoV (NC_004718) 
———————— BICoV-HKUS (NC_009021) 
4 0,55; 
1) MHV (NC_006852) 
0,96; PHEV (NC_007732) 
“HCoV-0C43 (NC_005147) 
1BtCoV (EF544563) 
BICoV (EF544565) 


1 4 ||BtCoV-HKU4 (0Q074652) 80,6-82,8 
a 78,3-78,6 


Pasa 


x $$ 4B ICovA 184049) 
0.99 VB (HQ184050) 


1; BtCoV (DQ648838) 


057 


PEDV (NC_003436) 
4;BtCoV (EU375862) 
BtCoV (EU375859) 
lexcov (€U375858) 
*BtCoV (EU375855) 
[ BtCov (EU375863) 
BtCoV(EU375861) 
| |BtCoV (EU375856) 
'Btcov (EU375857) 
CoV (EU375853) 
iT V (EU375854) 
| | ‘Btcov (EU375865) 
0,61 BtCoV (EU375860) 
rBtCoV (EU375869) 
BtCoV (EU375864) 
\BtCoV (EU375870) 
['BtCov (EU375868) 
077 4 BtCoV (EU375867) 
“BICoV K (HQ184060) 
BtCoV(DQ648822) 
BrCov (00648821) 
078 1'BxCov (00648824) 
al “BtCoV (DQ648823) 
(BtCoV (EU375875) 
!Btcov (EU375873) 


BACov (EU37S874) 
0,79 1| -BtCov (EU375872) 


{\BtCov (EU375866) 
0.8. .BICoVH (HQ184956) 
** | BtCoV (EU375871) 
0,98 BtCov (00648833) 
BtCoV-HKUS (DQ249224) 
BtCoV (DQ648837) 
4BtCoV-HKU2(0Q249235) 
BtCoV-HKU2 (DQ249213) 
rc BICoV (DQ648840) 
1 BICoV-HKUS (DQ249228) 
JBtCoV (EU834954) 
10.89 | BtCoV (EU834953) 
1 ‘BiCov (EU834952) 


0,59 0,92 


| BtCoV-HKU7 (DQ249226) 
—aA BICoV L (HQ184961) 
BICoV C (HQ184951) 
BICoV F (HQ184954) 
BICoV E (HQ184953) 
BCoV D (HQ184952) 
0,95, 0,92/BICoV G (HQ184955) eaieaa 
0.66 BICoV | (HQ184957) 
BICoV I" (HQT84958) 
Lt BtCoV (EU834951) + 801 
0,96 ——-HCoV-NL63 (NC_002645) 
HCoV-229E (NC_005831) 
1, FIPV (AY994055) 
—TGEV (NC_002306) 
IBV (NC001451) 


0.2 


HCoV-HKU1 (NC_006577) 


1/BtCoV (00648855) 
BICoV (00648854) 


Beta 


}o05 
} 73,2-74,3 


(~ 93,2-96,7 


rp 81,3-81,6 


} 94,7 Alpha 


f 95,8-96,6 


78,1-78,3 
74,5-76, 


75,8-76,1 


| Gamma 
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Fig. 2 CoV phylogenetic reconstruction based on 396 bp of the RdRp 
gene including 14 Spanish CoV from different bat species and 77 alpha, 
beta and gammacoronaviruses obtained from GenBank. Accession 
numbers are shown in brackets. BatCoV detected in Spain are 
highlighted in italics. For the analyses GTR substitution model, gamma 
estimation and two simultaneous runs of 10’ generations were done, 
each with four Markov chains, and the trees were sampled every 100 
generations. First 25% trees were excluded as burn-in from the analysis. 
Significant posterior probabilities are indicated. Complementary 
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information about sequences used in this phylogenetic reconstruction 
are shown in Table 2. Positive samples described in this work are shown 
in shaded rectangles and ovals. The new alpha and betacoronavirus 
groups described in this work are shown in shaded ovals. Amino acid 
identity was calculated with MEGA 4 using the pairwise deletion 
option. The alignment comprised the same 396 bp of the RdRp gene 
used for the phylogenetic reconstruction. Amino acid identities across 
132 amino acids are indicated next to the brackets that links every 
Spanish BatCoV to the phylogenetic neighbours 
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related to cluster HKU7 described in China [10]. BatCoV 
C, D, E, F, G and I belonged to the same lineage and 
showed an intriguing new independent cluster (significant 
posterior probability = 0.95) including BatCoV I’ (Fig. 2). 
In addition, BatCoV J and M were genetically related to [ 
CoV, although they did not really associate with any of 
their previously described linages. Sequences correspond- 
ing to two different groups of CoV were found in the same 
bat species (Hypsugo savii) as it was also found for other 
bat species elsewhere [12]. Spanish BatCoV described here 
were also classified calculating amino acid distances of 
these viruses from phylogenetic neighbours and related 
reference species [21]. Amino acid distance criteria 
recently described for separating RdRp grouping units 
(RGU) were adopted [21]. Interestingly, the amino acid 
distance criteria (>4.8% for alphacoronavirus and >6.3% 
for betacoronavirus) reinforced the presence of the new 
group of alphacoronaviruses mentioned above including 
Spanish BatCoV C, D, E, F, G and I; and additionally 
showed four new tentative groups (Fig. 2). BatCoV A and 
B were included in a new alphacoronavirus group and 
BatCoV I’ constituted another independent alphacorona- 
virus group. BatCoV J, and M represented two new beta- 
coronavirus groups and BatCoV H, K and L remained as 
part of several established species (Fig. 2). 

It has been previously suggested that some CoV asso- 
ciate to certain bat species [10, 11]. However, we found 
that different bat species from the same colony or location 
harbour CoV of the same genetic lineage (BatCoV A and 
B; G and J), indicating a greater diversity and higher 
complexity than previously described for the ecology of 
BatCoV. Similar exceptions were found in China and 
South America [12, 16, 26] and could also be observed 
with Australian BatCoV (Fig. 2 and Table 2). 

In conclusion, previous studies showed the presence of 
BatCoV in Europe. However, to our knowledge, this is the 
first report describing the presence of CoV RNA in Iberian 
bat species. Phylogenetic data indicate high diversity, wide 
distribution and complex ecology of CoV in bats captured 
in diverse Spanish locations. The sequences reported herein 
provide new insights into the diversity of coronaviruses and 
describe new phylogenetic linages that appear to diverge 
from all genotypes previously detected in other European 
locations. Future studies should clarify whether such 
apparently high diversity reflects the bio-geographical 
peculiarities of the Iberian Peninsula or not. This study 
contributes with a new dataset to the global surveillance of 
emerging BatCoV with pathogenic potential in humans. 
Our data reinforce the notion that the ecology and trans- 
mission of CoV in bat reservoirs is far from being com- 
pletely understood and that more studies will be necessary 
to evaluate the magnitude of the potential threat that these 
viruses pose to human health. 
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